Creep motion of a granular pile induced by thermal cycling 
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We report a time-resolved study of the dynamics associated with the slow compaction of a granular 
column submitted to thermal cycles. The column height displays a complex behavior: for a large 
amplitude of the temperature cycles, the granular column settles continuously, experiencing a small 
settling at each cycle; By contrast, for small-enough amplitude, the column exhibits a discontinuous 
and intermittent activity: successive collapses are separated by quiescent periods whose duration is 
exponentially distributed. We then discuss potential mechanisms which would account for both the 
compaction and the transition at finite amplitude. 



Literature canonically presents granular materials as a 
bunch of athermal particles P, Indeed, the energy 

necessary for the grains to hop one over another (a few 
tenths of millimeter, for instance) is roughly 10 orders 
of magnitude higher than the ambient thermal agitation 
ksT. This is probably one of the reasons why the behav- 
ior of a granular assembly submitted to temperature fluc- 
tuations has received so few attentions. Nonetheless, un- 
controlled thermal dilations of a granular pile have been 
reported to generate stress fluctuations large enough to 
hinder reproducible measurements of the stress fleld in- 
side the pile 0, 0| , and even suspected to be the driving 
factor leading to large-scale "static avalanches" Q. In- 
deed, the slow relaxation and compaction of a granular 
material, which has been hitherto produced by the input 
of mechanical energy fi\ , can be induced by periodically 
raising and then lowering the temperature of the gran- 
ular assembly, as recently brought to the fore by Chen 
and co-workers [sj. However, the compaction dynamics 
as well as the basic mechanisms at stake remain unknown. 
In this letter, we address the following questions: What 
is the dynamic of the top of a granular column submit- 
ted to thermal cycling? Does this compaction process 
exhibit features analogous to aging, as other compaction 
processes (tapping, cyclic shear...) do? And finally, what 
is the behaviour of the grain assembly in the limit of low 
amplitude temperature-cycles, i.e. well below a cycling 
amplitude of 40°C [1]? 

Compaction of granular materials has been thoroughly 
studied over the past 20 years both experimentally and 
theoretically It has been achieved by direct input of 
mechanical energy. By far, the most-used device consists 
of a vertical column which is periodically shaken with a 
controlled vibration intensity. Under vibration, a low- 
density granular-pile compacts towards a steady state 
with a final density that depends on the intensity of the 
vibrations [1, 2]. Nonetheless, those final states have been 
shown not to depend on the initial conditions (packing 
fraction) and even to be genuine thermodynamic states 
0. In addition to their general interest 0, vibrated 
granular materials have been found to exhibit deep com- 
mon physical features with out-of-equilibrium systems. 



such as glasses. First, the spectrum of density fluctua- 
tions was found to be strongly non-Lorentzian 10], which 
is a signature of multiple timescales in the system. Also, 
the response of the column to an abrupt change in shak- 
ing intensity is analogous to aging in glasses llj, . Fi- 
nally, in the limit of small-intensity vibrations, the relax- 
ation rate becomes so slow that the system cannot reach 
the steady-state density within the experiment timescale. 
At the microscopic scale, a grain exhibits random motion 
within confined volumes or cages, occasionally a longer 
exploration corresponding to a change of cage, and rarely 
jumps, during which it moves "sig nificantly" [13|. Cychc 
shear deformations of a granular assembly also lead to 
compaction and present the transient cage effects remind- 
ing one of colloidal glasses [3] ■ Those experimental 
results have motivated numerous theoretical and numer- 
ical works, a review of which can be found in [Ti] . 
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FIG. 1: Sketch of the experimental setup. Inset: pic- 
ture of the upper part of the column. The granular level is 
indicated by the white dotted-line. 

Up to now, the heating of a granular pile has only been 
used to develop dynamic light-scattering experiments [Tsj 
or to probe numerically its impact on the force network 
[l^ . Here, we report the time-resolved compaction dy- 
namics of a granular column submitted to a large number 
of temperature cycles and shed new lights on the key in- 
gredients governing this process. 
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Experimental set-up - The experimental set-up (Fig.[T]) 
consists in a vertical glass tube (height 1.7 m, inner diam- 
eter 13 mm) firmly fastened to a wall into the basement 
of the physics department, in order to flee mechanical vi- 
brations. The sample consists of a column (height H) of 
spherical glass beads (diameter d) poured into the tube. 
An additional gas- input, at the bottom of the column, 
makes it possible to loosen the pile thanks to an upward 
flow of dry nitrogen. The temperature cycles are im- 
posed by means of a heating cable (Prolabo, 40 W/m) 
directly taped on the outer surface of the tube wall. The 
resulting temperature is measured by means of a sen- 
sor (PtlOO, located close to the free surface of the gran- 
ular material) and a multimeter (Keithley Intruments, 
Model 196). The free surface of the material, which is 
illuminated by a red LED (Kingbright, L-793SRC-E, lo- 
cated inside the tube, above the granular material) is 
imaged from the side with a video camera (Panasonic, 
WV-BP500) connected to a frame grabber board (Data 
Translation, DT2255). A macro, running under a data- 
processing software (WaveMetrics, IGOR Pro 4.0), makes 
it possible to drive the heating power, to record the re- 
sulting variations of the temperature and to measure ac- 
curately the height H from the images: A subpixel reso- 
lution (namely, less than a tenth of a pixel which typically 
stands for 5 /xm) is achieved by considering the average 
position of the free surface, assumed to correspond to the 
inflection point in the vertical intensity-profile averaged 
over the whole diameter of the tube. Measurements are 
performed 20 times per temperature cycle. 

Due to long experimental times, we limit our report to 
a given diameter d = 510±90 /im of the grains (Matrasur 
Corp.) and to a given period 27r/ijj = 600 s of the cy- 
cles. The cycling period, 10 minutes, is arbitrarily chosen 
to be small enough to avoid excessively-long experimen- 
tal times but large enough to insure that the associated 
thermal penetration-length Ip = y^2X/ {p Cuj) ~ 6 mm 
is about the tube radius (A ~ 0.2 W m^^ K^^ and 
C ~ 10^ J m^"^ K^^ respectively denote the thermal 
conductivity and heat capacity of a typical glass-grains 
pile [13] ■) It is here crucial to note that the column is 
heated homogeneously along its whole length but that 
the temperature is likely to vary in the radial direction. 

Prior to each experiment, the granular column is pre- 
pared in a low-density state thanks to the dry-nitrogen 
upward flow. The top of the column is then higher than 
the field imaged by the camera (typically 1 cm above) 
and we set the amplitude of the cycles, AT, to the largest 
accessible value, AT = 27.1° C. The preparation of the 
sample ends when the top of the column enters the ob- 
servation field. At this point, the granular column is 
"quenched": The amplitude of the cycles is set to the 
chosen value AT lying bewteen and 27.1°C, which de- 
fines the origin of time t — i). The granular column is 
subsequently submitted to at least 1000 cycles (7 days). 

First experimental observations - Under the action of 
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FIG. 2: Height variation hn vs. number of cycles n. 

One observes first an exponential behavior at short time fol- 
lowed by a subsequent logarithmic creep at long time (The 
black curve corresponds to the test function defined in the 
text.) Inset: Oscillations of the column height associated with 
the temperature cycles [An and (5„ are respectively defined to 
be the amplitude of the increase and the drift of /i„ at the 
cycle n.) For the chosen AT, the column settles slightly at 
each cycle {H = 140 cm, 2-k/uj = 600 s and AT = 10.8° C.) 



the temperature cycles the column height decreases: we 
report the variation /i„ = H{2iTn/Lo) — H{0), where n de- 
notes the time in units of the cycle period or, equivalently 
when integer, the number of imposed cycles (Fig. [2]). We 
observe that the thermal-induced compaction is a very 
slow phenomenon: after 7 days (1000 cycles), the de- 
crease of the height is of about 1.5 cm (about 1% of the 
height H), which indicates that the system remains very 
far from the maximum compaction (roughly a decrease 
of about 10% of the column height, i.e. /i„ ~ 10 cm) 
within the experimental time. Accordingly, after the ex- 
periment, we checked that a single finger tap produces 
a collapse of the granular height of one centimeter at a 
rough estimate. We also checked during 3 days that the 
height H of the column stays constant when no temper- 
ature cycles are imposed, which proves that ambient me- 
chanical vibrations and changes in the room temperature 
have no (or little) effect in our experimental conditions. 

Before we discuss the compaction dynamics, it is rel- 
evant to determine which part of the column is involved 
in the phenomenon. Our measurements are accurate 
enough for observing oscillations of H associated with 
the temperature variations (Fig. ^ inset). We observe 
on the raw data that the amplitude An , which is propor- 
tional to AT, increases logarithmically with n (Fig. [3]). 
The oscillations of H are due to the thermal dilation of 
both the tube and the granular material. In order to 
assess the contribution of the granular material, we first 
determine the amplitude, 6ht{z), of the tube displace- 
ment in the laboratory frame as a function of the height 
z (origin at the bottom of the column) by marking its 
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FIG. 3: Amplitude An vs. number of cycles n. The data 
are successfully accounted for by A„ = AT [ao + foo lii(n)] with 
ao ~ 13.8 nm/K and bo ~ 2.5 nm/K {H = 140 cm, 2tv/uj = 
600 s and, from bottom to top, AT = 10.8, 16.2 and 27.1° C.) 



outer wall. The amplitude 5ht {z) is found to be linear in 
z and the slope provides us with an estimate of the linear 
thermal-expansion coefficient, k = (3.6±0.4) xlO^^ K"-'^, 
of the tube material. Then, considering the relative vari- 
ation of the inner volume, we write the relation between 
the amplitude A„ and the relative variation, SVg/Vg, of 
the volume Vg of the granular material: An — 5ht{0) = 
H{5Vg /Vg — 2kAT). Experiments performed for different 
height H demonstrate that A^ — Sht(0) is proportional 
to H, which shows that SVg/Vg is independent of H. We 
thus infer that the whole height H of the granular column 
is involved in the observed oscillations of the free surface. 

Compaction dynamics - For a large AT (typically more 
than 3°C), the column systematically compacts during 
the decrease of the temperature, at each cycle, and the 
amplitude (5„ of the collapses (defined in Fig. [2l inset) 
decreases when n increases, which accounts for the over- 
all behavior of as a function of n. We show here 
that the results obtained by Chen and co-workers for 
AT > 40°C [i| hold true at much smaller AT. Besides, 
we resolve the column dynamics and, this, over a sig- 
nificant longer duration. Here, the column height first 
decreases exponentially to reach, after the characteris- 
tic number Uc of cycles, a logarithmic behavior at long 
time (Fig. [5]). This response to the thermal quenching 
is very similar to the one the system exhibits to step 
strain solicitation . We estimated ric from the exper- 
imental data by interpolating ft,„ with the test function 
/i^ = Hq + he exp(— n/ric) -I- hi ln(n). In addition, we 
checked that the shortest characteristic time obtained 
from the interpolation by a double exponential [31 is 
consistent with ric and presents the same behavior as a 
function of AT. We observe that ric drastically increases 
when AT is decreased. Moreover, nc even seems to di- 
verge at ATc (We get 2.8 < AT^ < 4.0 °C, Fig.gD, which 
prefigures a qualitatively-different dynamics of the col- 
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FIG. 4: Characteristic number ric vs. amplitude 

AT. The characteristic number of cycles increases dras- 
tically when AT is decreased and even seems to diverge for 
ATc ~ 3°C. Inset: For AT < AT^, the column settles, 
by jumps, linearly with time. For readibility, we only display 
data obtained during the first half of the experiment duration 
(14 days, H = 140 cm, 2n/io = 600 s and AT = 2.8°C.) 



umn compaction at a lower amplitude. 

For AT < ATc, one indeed observes that the column is 
not flowing regularly anymore, but evolves by successive 
collapses separated by rest periods (Fig. IH inset). The 
overall compaction velocity is constant (steady regime) 
and it is therefore relevant to consider the probability 
distribution P{An) of the number of cycles An between 
two successive collapses (Fig. [S]). However, due to the 
small number of collapses, TV = 175, within the experi- 
mental time (14 days), the behavior of P(An) is difficult 
to assess directly and we consider the cumulated prob- 
abihty, F{An) = J^^Pix) dx, instead [23|. We obtain 
F{An) ~ exp(-An/An*) with An* ~ 12 for AT = 
2.8°C (Fig. [S^). The delay An between two successive 
collapses is thus observed to be exponentially distributed. 
By contrast, the amplitude of the corresponding col- 
lapses, Sn, exhibits a Gaussian probability-distribution 
Q{S) whose width is estimated to be S* = 63 ± 4/im 
(Fig. Oj), roughly a tenth of the grain diameter. The 
compaction of the material is thus shown to result from 
collapses that are randomly distributed in time. 

Discussion and conclusion - The mechanisms leading 
to the compaction could be either extrinsic (the dilation 
of the container at stake) or intrinsic (only due to the 
dilation of the grains.) Chen et al. suggested that the 
compaction could be due to "the difference between the 
thermal expansion of the container and of the grains" @] , 
thus to an extrinsic mechanism. In our case, the beads 
(Kf, = 9 X 10-6 K-i) and the container (k = 3.6 x 10"*^ 
K^^) are made of glass and present roughly the same 
thermal expansion coefficient 2l|. Thus, we do not think 
that the thermal dilation of the container is the primary 
cause of the compaction. This conclusion is supported 
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FIG. 5: Statistics below the transition, (a) : The cumu- 
lated probability F{An) ~ exp (—An/An*) is an exponen- 
tial function of An and the characteristic number An* ~ 12. 
(b) : The probability distribution Q{5) ~ exp [— (5/5*)^] is 
a Gaussian function whose width S* = 63 ± 4 ^m. The sta- 
tistical data reported herein correspond to the 175 collapses 
observed within the experimental time (14 days, H — 140 cm, 
27r/a; = 600 s and AT = 2.8°C). 



by preliminary results obtained by heating periodically 
the granular column along its axis 12211 and the recent 
simulations by Vargas and McCarthy |16| . Dealing with 
an intrinsic mechanism, we guess that an overall tem- 
perature change would only produce a homogeneous di- 
lation of the pile and no spatial reorganisation of the 
grains. We thus propose the following scenario, compat- 
ible with previous results due to the radial thermal- 
gradient of amplitude AT/lp, the grains are subjected to 
a periodic shear, 7, along the vertical which induces lo- 
cal rearrangments and, accidentally, internal avalanches. 
This mechanism is supported by the estimate of the shear 
7 = Ki,{AT/lp)L, where L is a characteristic size which 
ranges from the grain- to the tube-diameter. We obtain, 
for L £ [0.510, 13] mm, that the relative displacement 
of two neighbor grains F = ^d/AT e [0.4,10] nm/°C 
whereas, using an AFM, we estimated the maximum size 
of the asperities associated with the surface roughness 
of the beads, s ~ 100 nm (Invoking the surface rough- 
ness is supported by numerical simulations: for perfectly 
smooth grains the maximum compacity is reached after 
a few cycles [l^, which contrasts with our observations.) 
In this framework, the continuous compaction would thus 
correspond to amplitudes AT such as FAT > s, whereas 
jumps would be observed for AT < ATc ~ s/F. From 
the experimental ATc — 3°C, we could infer that the 
typical length at stake, L, is about the tube radius and, 
thus, the transition a finite-size effect. Due to long ex- 
perimental times, we limited our report to results which 
clearly demonstrate the efficiency of the techniques. The 
dependance of the critical amplitude ATj, on the column 



height and diameter as well as on the cycling frequency 
and grain roughness is necessary to corroborate or dis- 
card the proposed mechanism. In particular, the sensi- 
tivity to the frequency uj will indicate the role played by 
the temperature gradient. We also plan to locate the ori- 
gin of the collapses by listening to the associated noise 
and to observe the local displacement by optical means. 
This way, we hope to correlate the macroscopic behavior 
to the local dynamics as already done for tapping l^, l^] 
and cycling shear experiments [3, IIS] • 
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